A magnetic origin for high temperature superconductivity in iron pnictides 
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In conventional Bardeen-Cooper-Schrieffer (BCS) superconductors 1 1 1, superconductivity occurs when electrons form 
coherent Cooper pairs below the superconducting transition temperature T^. Although the kinetic energy of paired elec- 
trons increases in the superconducting state relative to the normal state, the reduction in the ion lattice energy is sufficient 
to give the superconducting condensation energy (Ec = -N(0)A^/2 and A ^ 2^^^)/)^"^^^^^^^% where N(0) is the electron den- 
sity of states at zero temperature, Hcod is the Debye energy, and Vq is the strength electron-lattice coupling) | 2 , 3 , 13 1. For 
iron pnictide superconductors derived from electron or hole doping of their antiferromagnetic (AF) parent compounds 
dHIOl, the microscopic origin for superconductivity is unclear llIT! . Here we use neutron scattering to show that high-r^ 
superconductivity only occurs for iron pnictides with low-energy (< 25 meV or ~ 6.5^^7^) itinerant electron-spin exci- 
tation coupling and high energy (> 100 meV) spin excitations. Since our absolute spin susceptibility measurements for 
optimally hole-doped iron pnictide reveal that the change in magnetic exchange energy below and above Tc can account 
for the superconducting condensation energy, we conclude that the presence of both high-energy spin excitations giving 
rise to a large magnetic exchange coupling / and low-energy spin excitations coupled to the itinerant electrons is essential 
for high-r^ superconductivity in iron pnictides. 

For BCS superconductors, the superconducting condensation energy Ec and Tc are controlled by the strength of the Debye 
energy Hojd and electron-lattice coupling Vq | 1 |. A material with large Hojd and lattice exchange coupling is a necessary but 
not a sufficient condition to have high-r^ superconductivity. On the other hand, a soft metal with small Hojo (such as lead and 
mercury) will also not exhibit superconductivity with high-rc. For unconventional superconductors such as iron pnictides, the 
superconducting phase is derived from hole and electron doping from their AF parent compounds ItSHTOl . Although the static 
long-range AF order is gradually suppressed when electrons or holes are doped into the iron pnictide parent compound such 
as BaFe2As2 |6-T0'|, short-range spin excitations remain throughout the superconducting phase and are coupled directly with 
the occurrence of superconductivity |[4H6l [T2l [T3l [15- 191 . For spin excitations mediated superconductors, the superconducting 
condensation energy should be accounted for by the change in magnetic exchange energy between the normal (N) and super- 
conducting (S) phases at zero temperature. Within the t-J model |22|, AEex(T) = 2J[{Si+x • S/)^^ - (8/+^ • S/)^], where (8/+^ • S/) 
is the dynamic spin susceptibility in absolute units at temperature T jSWTSl. 

To determine how high-r^ superconductivity in iron pnictides is associated with spin excitations, we consider the phase 
diagram of electron and hole-doped iron pnictide BaFe2As2 (Fig. la) \TU\. In the undoped state, BaFe2As2 forms a metallic low- 
temperature orthorhombic phase with collinear AF structure as shown in the inset of Fig. la. Inelastic neutron scattering (INS) 
measurements have mapped out spin waves throughout the Brillouin zone, and determined the effective magnetic exchange cou- 
plings (2]. Upon doping electrons to BaFe2As2 by partially replacing Fe with Ni to induce superconductivity in BaFe2-xNi;cAs2 
with maximum ^ 20 K at Xg = 0.1 I24 L the low-energy (< 80 meV) spin waves in the parent compounds are broadened and 
form a neutron spin resonance coupled to superconductivity |[T6HT9l , while high-energy spin excitations remain unchanged O. 
With further electron-doping to Xe > 0.25, superconductivity is suppressed and the system becomes a paramagnetic metal (Fig. 
la) |24|. For hole-doped Bai_;cKjcFe2As2 li25J , superconductivity with maximum Tc = 38.5 K appears at Xh ~ 0.33 |6| and pure 
KFe2As2 at = 1 is a = 3.1 K superconductor Q. In order to determine how spin excitations throughout the Brillouin zone 
are correlated with superconductivity in iron pnictides, we study optimally hole-doped Bao.67Ko.33Fe2As2 (Tc = 38.5 K, Fig. 
Ic), pure KFe2As2 (Tc = 3 K, Fig. lb), and nonsuperconducting electron-overdoped BaFe1.7Nio.3As2 (Fig. Id). If spin excita- 
tions are responsible for mediating electron pairing and superconductivity, the change in magnetic exchange energy between the 
normal (N) and superconducting (S ) state should be large enough to account for the superconducting condensation energy fTT| . 

In previous INS work on powder |[T2l[T3l and single crystals |4 | of hole-doped Bai_;cKxFe2As2, low-energy spin excitations 
were found to be dominated by a resonance coupled to superconductivity. In fact, density functional theory (DFT) calculations 
based on sign reversed quasiparticle excitations between hole Fermi surface pocket near the Brillouin zone center and electron 
pocket near the zone comer (see supplementary information) 1 11 1 predict correctly the longitudinally elongated spin excitations 
from Qaf = (1,0) for optimally hole-doped iron pnictide (inset in Fig. la) El [191 . Foi" pure KFe2As2 (xh = 1), low-energy 
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(< 14 meV) spin excitations become longitudinally incommensurate from Qaf (inset in Fig. la) ifTSl . These results, as well as 
work on electron-doped iron pnictides that reveal transversely elongated spin excitations from Qaf OEIIISl, have shown that 
the low energy nematic-like spin dispersion in iron-based superconductors can be accounted for by itinerant electrons on hole 
and electron nested Fermi surfaces 1 10 |. 

Here we use INS to show that the effective magnetic exchange couplings in hole-doped Bao.67Ko.33Fe2As2 only soften slightly 
from that of AF BaFe2As2 EI and electron-doped BaFe1.9Nio.1As2 superconductor (Fig. If) (61. The effect of hole-doping 
in BaFe2As2 is to suppress high-energy spin excitations and transfer the spectral weight to low-energies that couple to the 
appearance of superconductivity (Fig. Ih). This is qualitatively consistent with theoretical methods based on DFT and dynamic 
mean filed theory (DMFT, see supplementary information) |7|. By using the INS measured magnetic exchange couplings and 
spin susceptibility in absolute units, we calculate the superconductivity-induced lowering of magnetic exchange energy and find 
it to be about seven times larger than the superconducting condensation energy determined from specific heat measurements 
for Bao.67Ko.33Fe2As2 ifTSl . These results are consistent with spin excitations mediated electron pairing mechanism |[T3l . For 
the nonsuperconducting electron overdoped BaFe1.7Nio.3As2, we find that while the effective magnetic exchange couplings only 
reduce slightly compared with that of optimally electron-doped BaFe1.9Nio.1As2 (Fig. Ig) |6|, the low-energy spin excitations 
(< 50 meV) associated with hole and electron pocket Fermi surface nesting disappear, thus revealing the importance of Fermi 
surface nesting and itinerant electron- spin excitation coupling to the occurrence of superconductivity (Fig. Ih). Finally, for 
heavily hole-doped KFe2As2 with low-J^ superconductivity (Fig. lb), there are only incommensurate spin excitations below 
~25 meV and the correlated high-energy spin excitations prevalent in electron-doped and optimally hole-doped iron pnictides 
are completely suppressed (Fig. le), indicating a dramatic softening of effective magnetic exchange coupling (inset Fig. Ih). 
Therefore, high-r^ superconductivity in iron pnictides requires two fundamental ingredients: a large effective magnetic exchange 
coupling ITSl, much like large Debye energy for high-r^ BCS superconductors, and a strong itinerant electrons- spin excitations 
coupling from Fermi surface nesting [11 1, like electron-phonon coupling in BCS superconductors. The presence of correlated 
electronic states exhibiting both local and itinerant properties is essential for the mechanism of superconductivity 1 10 |. 

To substantiate the key conclusions of Fig. 1, we present the two-dimensional (2D) constant-energy images of spin excitations 
in the {H, K) plane at different energies for KFe2As2 (Figs. 2a-2c), Bao.67Ko.33Fe2As2 (Figs. 2d-2f), and BaFe1.7Nio.3As2 (Figs. 
2g-2i) above T^. In previous INS work on KFe2As2, longitudinal incommensurate spin excitations were found by triple axis 
spectrometer measurements for energies from 3 to 14 meV in the normal state ifTSl . While we confirmed the earlier work using 
time-of-flight INS for energies below E = 15 + 1 meV (Figs. 2a and 2b), our new data collected at higher excitation energies 
reveals that incommensurate spin excitations converge into a broad spin excitation near E = 2{) meV and disappear for energies 
above 25 meV (Fig. 2c). Therefore, there are no measurable correlated spin excitations for energies above 30 meV, indicating 
that the effective magnetic coupling in KFe2As2 has reduced about 90% compared with that of BaFe2As2 (see supplementary 
information). For Bao.67Ko.33Fe2As2, spin excitations at £" = 5 ± 1 meV are longitudinally elongated from Qaf as expected from 
DFT calculations (Fig. 2d) ||4l[T9l. At the resonance energy (£" = 15 + 1 meV) [12J, spin excitations are isotropic above Tc (Fig. 
2e). On increasing energy further to £" = 50 + 10 meV, spin excitations change to transversely elongated from Qaf similar to 
spin excitations in optimally electron-doped superconductor BaFe1.9Nio.1As (Fig. 2f) |6|. Figures 2g-2i summarize similar 2D 
constant-energy images of spin excitations for nonsuperconducting BaFe1.7Nio.3As. At £" = 9 + 3 (Fig. 2g) and 30+10 meV 
(Fig. 2h), there are no correlated spin excitations near the Qaf- Upon increasing energy to £" = 59 + 10 meV (Fig. 2i), we see 
clear spin excitations transversely elongated from Qaf (Fig. 2i). 

Figures 3a-3d show 2D images of spin excitations in BaFe1.7Nio.3As2 at £^ = 70 + 10, 112 + 10, 157 + 10, and 214 + 10 meV, 
respectively. Figures 3e-3h show wave vector dependence of spin excitations at energies £^ = 70 + 10, 115 + 10, 155 + 10, and 
195 + 10 meV, respectively, for Bao.67Ko.33Fe2As2. Similar to spin waves in BaFe2As2 |2|, spin excitations in BaFe1.7Nio.3As2 
and Bao.67Ko.33Fe2As2 split along the ^-direction for energies above 80 meV and form rings around 2 = (+1, +1) positions near 
the zone boundary, albeit at slightly different energies. Comparing spin excitations in Figs. 3a-3d for BaFei.7Nio.3 As2 with those 
in Figs. 3e-3h for Bao.67Ko.33Fe2As2 in absolute units, we see that spin excitations in BaFe1.7Nio.3As2 extend to slightly higher 
energies and have larger intensity above 100 meV. 

To understand the wave vector dependence of spin excitations in hole-doped Bao.67Ko.33Fe2As2, we have carried out the 
random phase approximation (RPA) calculation of the dynamic susceptibility in a pure itinerant electron picture using method 
described before Q. Figures 3i and 3j show RPA calculations of spin excitations at £ = 70 and 155 meV, respectively, for 
Bao.67Ko.33Fe2As2 assuming that hole doping induces a rigid band shift |5 |. The outcome is in clear disagreement with Figs. 
3e and 3g, indicating that a pure RPA type itinerant model cannot describe the wave vector dependence of spin excitations in 
hole-doped iron pinctides at high energies. For comparison with the RPA calculation, we also used a combined DFT and DMFT 
approach to calculate the imaginary part of the dynamic susceptibility ;^'XQ. ^) in the paramagnetic state. Figures 3k and 
3i show calculated spin excitations at £ = 70 and 155 meV, respectively. Although the model still does not agree in detail with 
the data in Figs. 3e and 3g, it captures the trend of spectral weight transfer away from Qaf = (1, 0) on increasing the energy and 
forming a pocket centered at 2 = (1, 1). 

By carrying out cuts through the 2D images similar to Figs. 2d-2f and 3e-3h along the [1, ^] and [//, 0] directions (see supple- 
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mentary information), we establish the spin excitation dispersion along the two high symmetry directions for Bao.67Ko.33Fe2As2 
and compare with the dispersion of BaFe2As2 (Fig. If) L2J. In contrast to the dispersion of electron-doped BaFe1.9Nio.1As2 O, 
we find clear softening of the zone boundary spin excitations in hole-doped Bao.67Ko.33Fe2As2 from spin waves in BaFe2As2 El- 
We estimate that the eff'ective magnetic exchange coupling in Bao.67Ko.33Fe2As2 is reduced by about 10% from that of BaFe2As2 
(see supplementary information). Similarly, Figure Ig shows the dispersion curve of BaFeijNioj As2 along the [1, K] direction 
plotted together with that of BaFe2As2 EI- For energies below ~50 meV, spin excitations are completely gapped marked in the 
dashed area probably due to the missing hole-electron Fermi pocket quasiparticle excitations |1JJ1281. Based on the 2D spin 
excitation images similar to Figs. 2a-2c, we plot in Fig. le the dispersion of incommensurate spin excitations in KFe2As2. The 
incommensurability of spin excitations is weakly energy dependent below £" = 12 meV but becomes smaller with increasing 
energy above 12 meV (see supplementary information). Correlated spin excitations for energies above 25 meV are completely 
suppressed as shown in the shaded area in Fig. le. 

To quantitatively determine the eff'ect of electron and hole doping on the overall spin excitations spectra, we calculate the local 
dynamic susceptibility per formula unit (f.u.) in absolute units, defined as;^'^^) = / Af''(Q'^)^Q/ / where ;^''(q. ^) = 
(l/3)rr(Y^^(q, 61;)), at diff'erent energies for Bao.67Ko.33Fe2As2, BaFe1.7Nio.3As2, and KFe2As2. Figure Ih shows the outcome 
together with previous data on optimally electron-doped superconductor BaFe1.9Nio.1As2 O. While electron-doping up to 
BaFe1.7Nio.3As2 does not change the spectral weight of high-energy spin excitations from that of BaFe1.9Nio.1As2, hole-doping 
dramatically suppresses the high-energy spin excitations and shift the spectral weight to lower energies (Fig. Ih). For heavily 
hole-doped KFe2As2, spin excitations are confined to energies below about E = 25 meV (inset in Fig. Ih). The reduction 
of the high energy spin spectral weight and its transfer to low energy with hole doping, but not with electron doping, is not 
naturally explained by the band theory, and requires models which incorporate both the itinerant quasiparticles and the local 
moment physics (see supplementary information) |7 |. The hole doping makes electronic state more correlated, as local moment 
formation is strongest in the half-filled shell, and mass enhancement larger thereby reducing the electronic energy scale in the 
problem. 

Finally, to determine how low-energy spin excitations are coupled to superconductivity in Bao.67Ko.33Fe2As2, we carried out a 
detailed temperature dependent study of spin excitation atE = 15+1 meV. Comparing with strongly c-axis modulated low-energy 
(E <7 meV) spin excitations [4 ], spin excitations at the resonance energy are essentially 2D without much c-axis modulations. 
Figures 4a-4d show the 2D mapping of the resonance at T = 25, 38,40, and 45 K, respectively. While the resonance reveals 
a clear oval shape at temperatures below Tc (Figs. 4a and 4b), it changes into an isotropic circular shape abruptly at Tc (Figs. 
4c and 4d) as shown by the dashed lines representing full-width-at-half-maximum (FWHM) of the excitations. Temperature 
dependence of the resonance width along the [//, 0] and [1,^] directions in Fig. 4e reveals that the isotropic to anisotropic 
transition in momentum space occurs at Tc. Figure 4f shows temperature dependence of the resonance from 9 K to 40 K, 
which vanishes at Tc. Figure 4g plots temperature dependence of the mode energy together with the sum of the superconducting 
gaps from hole and electron pockets 1281 . Figure 4h compares temperature dependence of the superconducting condensation 
energy 1 15 1 with superconductivity induced intensity gain of the resonance. By calculating spin excitations induced changes in 
magnetic exchange energy (see supplementary information) 1 13], we find that the diff'erence of magnetic exchange interaction 
energy between the superconducting and normal state is approximately seven times larger than the superconducting condensation 
energy 1 15|, thus identifying AF spin excitations as the major driving force for superconductivity in Bao.67Ko.33Fe2As2. 

One way to quantitatively estimate the impact of hole/electron doping and superconductivity to spin waves of BaFe2As2 is to 
determine the energy dependence of the local moment and total fluctuating moments (m^^ IH. From Fig. Ih, we see that hole- 
doping suppresses high-energy spin waves of BaFe2As2 and pushes the spectral weight to resonance at lower energies. The total 
fluctuating moment of Bao.67Ko.33Fe2As2 below 300 meV is (rrp-^ = l.l ± 0.3 per Fe, somewhat smaller than 3.2 ± 0. 16 yu^ per Fe 

forBaFe2As2 andBaFe1.9Nio.1As2 O. For comparison, BaFe1.7Nio.3As2 andKFe2As2 have (m^^ = 2.74 + 0.11 and 0.1 +0.02 
per Fe, respectively. Therefore, the total magnetic spectral weights for difl'erent iron pnictides have no direct correlation with 
their superconducting T^'s. From Fig. Ih, we also see that the spectral weight of the resonance and low-energy (< 100 meV) 
magnetic scattering in Bao.67Ko.33Fe2As2 is much larger than that of electron-doped BaFe1.9Nio.1As2. This is consistent with 
a large superconducting condensation energy in Bao.67Ko.33Fe2As2 since its efl'ective magnetic exchange coupling / is only 
-10% smaller than that of BaFe1.9Nio.1As2 (Fig. Ih) |[T5][29l. For electron-overdoped nonsuperconducting BaFe1.7Nio.3As2, the 
lack of superconductivity is due to the absence of low-energy spin excitations coupled to the hole and electron Fermi surface 
nesting even though the efl'ective magnetic exchange remains large IHlim. This means that by eliminating [(8/+^ • S/)^ - 
{Si+x • S/)^ ], there is no magnetic driven superconducting condensation energy and thus no superconductivity. On the other hand, 
although the complete suppression of correlated high-energy spin excitations in KFe2As2 can dramatically reduce the efl'ective 
magnetic exchange coupling in KFe2As2 (Fig. le), one can still have superconductivity with reduced Tc. If spin excitations are 
a common thread of the electron pairing interactions for unconventional superconductors 1 13 1, our results reveal that both the 
large efl'ective magnetic exchange couplings and itinerant electron-spin excitation interactions are essential ingredients to achieve 
high-Tc superconductivity, much like large Debye energy and strength of electron-lattice coupling are necessary for high-r^ 
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BCS superconductors. Therefore, the presence of correlated electronic states exhibiting wave-particle duality is important for 
the mechanism of high-r^ superconductivity in iron pnictides. 
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FIG. 1: Summary of transport, magnetic, and neutron scattering results. Our experiments were carried out on the MERLIN and 
MAPS time-of-flight chopper spectrometers at the Rutherford- Appleton Laboratory, UK |i2n6J. We co-aUgned 19 g of single crystals of 
Bao.67Ko.33Fe2As2 (with in-plane and out-of-plane mosaic of ~4°), 3 g of KFe2As2 (with in-plane and out-of-plane mosaic of ~7.5°), and 40 
g of BaFe1.7Nio.3As2 (with in-plane and out-of-plane mosaic of ~3°). Various incident beam energies were used as specified, and mostly with 
Ei parallel to the c-axis. To facilitate easy comparison with spin waves in BaFe2As2 |2|, we defined the wave vector Q at {q^, qy, qz) as 
(H,K,L) = (qxa/27T,qyb/2n,qzC/27T) reciprocal lattice units (rlu) using the orthorhombic unit cell, where a - b - 5.51 A, and c - 13.135 
A for Bao.67Ko.33Fe2As2, a = b = 5.43 A, and c = 13.8 A for KFe2As2, and a = Z? = 5.6 A, and c = 12.96 A for BaFe1.7Nio.3As2. The 
data are normalized to absolute units using a vanadium standard with 20% errors |6|. (a) The electronic phase diagram of electron and 
hole-doped BaFe2As2 1 10|. The right inset shows crystal and AF spin structures of BaFe2As2 with marked the nearest (Jia, Jib) and next 
nearest neighbor (^2) magnetic exchange couplings. The left insets show the evolution of low-energy spin excitations in Bai_;cK;(:Fe2As2. (b,c) 
Temperature dependence of magnetic susceptibility for our KFe2As2 and Bao.67Ko.33Fe2As2. (d) Temperature dependence of the resistivity for 
BaFe1.7Nio.3As2. (e,f,g) The filled circles are spin excitation dispersions of KFe2As2 at 5 K, Bao.67Ko.33Fe2As2 at 9 K, and BaFe1.7Nio.3As2 
at 5 K, respectively. The shaded areas indicate vanishing spin excitations and the solid lines show spin wave dispersions of BaFe2As2 lU. 
(h) Energy dependence ofx'\(^) for BaFe1.9Nio.1As2 (dashed line), BaFe1.7Nio.3As2 (green solid circles), Bao.67Ko.33Fe2As2 below (Solid red 
circles and solid red line) and above (open purple circles and solid lines) T^. The inset shows Energy dependence of ;^''(^) for KFe2As2. The 
vertical error bars indicate the statistical errors of one standard deviation. The horizontal error bars in (h) indicate the energy integration range. 
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FIG. 2: Constant-energy slices through magnetic excitations of KFe2As2, Bao 67K0 33Fe2As2, and BaFei vNio 3AS2 at different energies. 

The color bars represent the vanadium normahzed absolute spin excitation intensity in the units of mbarn/sr/meV/f.u. Two dimensional images 
of spin excitations at 5 K for KFe2As2 (a) £" = 8 + 3 meV obtained with Ei - 20 meV. The right side incommensurate peak is obscured by 
background scattering, (b) 13 + 3 meV with Ei - 35 meV, (c) 53 + 10 meV with Ei - 80 meV. For Bao.67Ko.33Fe2As2 at T = 45 K, images of 
spin excitations at (d) £^ = 5 + 1 meV obtained with Ei - 20 meV, (e) 15 + 1 meV with Ei - 35 meV, and (f) 50 + 2 meV obtained with Ei - 80 
meV. The dashed box in (d) indicates the AF zone boundaries for a single FeAs layer and the black dashed lines mark the orientations of spin 
excitations at different energies. Images of spin excitations for BaFe1.7Nio.3As2 at T = 5 K and {g) E - 9 ±3 meV obtained with Ei - 80 meV, 
(h) 30+10 meV with Ei - 450 meV, and (i) 59 + 10 meV with Ei - 250 meV. The white crosses indicate the position of Qaf- 
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FIG. 3 : Constant-energy images of spin excitations in the 2D [H, K] plane at different energies for BaFei vNio 3 As2 and Bao 67K0 33Fe2 As2 
and its comparison with RPA/DMFT calculations for Bao.67Ko.33Fe2As2. Spin excitations of BaFe1.7Nio.3As2 at energy transfers of (a) 
£ = 70 + 10 meV obtained with Ef = 250 meV;(b) 112+10 meV Et = 250 meV;(c) 157 + 10 meV; (d) 214 + 10 meV, all obtained with 
Ei = 450 meV at 5 K. A flat backgrounds have been subtracted from the images. Spin excitations of Bao.67Ko.33Fe2As2 at energy transfers of 
(e) ^ = 70 + 10 meV obtained with Et = 170 meV;(f) 115 + 10 meV;(g) 155 + 10 meV; (h) 195 + 10 meV, all obtained with Et = 450 meV at 9 
K. Wave vector dependent backgrounds have been subtracted from the images. RPA calculations | 5 1 of spin excitations for Bao.67Ko.33Fe2As2 
at (i) E = 10 meV and G) ^ = 155 meV. DMFT calculations for Bao.67Ko.33Fe2As2 at (k) E = 10 meV and (1) E = 155 meV. 
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FIG. 4: Wave vector, energy, and temperature dependence of the resonance across Tc = 38.5 K. Constant-energy (E = 15 ± \ meV) 
images of spin excitations at (a) T = 25, (b) 38, (c) 40, and (d)45 K obtained with Ei = 35 meV. In order to make fair comparison of the 
scattering line shape at different temperatures, the peak intensity at each temperature is normaUzed to 1. The pink and green arrows in (a) 
mark wave vector cut directions across the resonance. The integration ranges are -0.2 < ^ < 0.2 along the [H, 0] direction and 0.8 <H< 1.2 
along the [1, ^] direction. The FWHM of spin excitations are marked as dashed lines, (e) The FWHM of the resonance along the [H, 0] and 
[1, ^] directions as a function of temperature across Tc. (f) Energy dependence of the resonance obtained by subtracting the low-temperature 
data from the 45 K data, and correcting for the Bose population factor, (g) The black diamonds show temperature dependence of the sum of 
hole and electron pocket electronic gaps obtained from Angle Resolved Photoemission experiments for Bao.67Ko.33Fe2As2 | 28 1. The red solid 
circles show temperature dependence of the resonance, (h) Temperature dependence of the superconducting condensation energy from heat 
capacity measurements |15| and the intensity of the resonance integrated from 14 to 16 meV. The error bars indicate the statistical errors of 
one standard deviation. 
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Supplementary Information: A magnetic origin for high-temperature superconductivity in iron pnictides 

Meng Wang*, Chenglin Zhang*, Xingye Lu*, Guotai Tan, Huiqian Luo, Yu Song, Miaoyin Wang, Xiaotian Zhang, E. A. Gore- 
mychkin, T. G. Perring, T. A. Maier, Zhiping Yin, Kristjan Haule, Gabriel KotHar, Pengcheng Dai 

ADDITIONAL DATA AND ANALYSIS 

We first discuss detailed experimental results on electron-doped iron pnictides, focusing on comparison of electron over- 
doped nonsuperconducting BaFe1.3Nio.3As2 with optimally electron-doped superconductor BaFe1.9Nio.1As2 and antiferromag- 
netic (AF) BaFe2As2. SFigure[5]shows the evolution of Fermi surfaces as a function of increasing Ni-doping obtained from the 
tight-binding model of Graser et al. 1 1 1 and our measured low-energy spin excitations spectra for BaFe2As2, BaFe1.9Nio.1As2, 
and BaFe1.3Nio.3As2. The absence of hole Fermi pocket near the zone center for BaFe1.3Nio.3As2 means quasiparticle excitations 
between the hole and electron pockets are not possible, thus eliminating low-energy spin excitations. SFigure[6]shows detailed 
comparison of spin excitations at diff'erent energies for BaFe1.9Nio.1As2 and BaFe1.3Nio.3As2. SFigure[7]plots the evolution of 
spin excitations for BaFe1.3Nio.3As2 in several Brillouin zones. SFigure [8] shows the identical wave vector versus energy cuts 
for BaFe2As2 | 2 | and BaFe1.7Nio.3As2. The spin excitations in BaFe1.7Nio.3As2 are clearly absent below 50 meV, which is much 
bigger than the 15 meV single ion spin anisotropy gap in BaFe2As2 El. SFigure [9] shows the comparison of cuts along two high- 
symmetry directions for BaFe2-xNi;cAs2 at Xe = 0, 0.1, 0.3. SFigure [T0| plots the constant-Q cuts and spin correlation lengths for 
these three samples. It is clear that the zone boundary spin excitations for electron-doped materials remain almost unchanged at 
least up to Xe = 03. 
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SFig 5: Schematics of Fermi surface evolution as a function of Ni-doping for BaFe2-jcNi;cAs2 and corresponding spin excitations. (a,b,c) 
Evolution of Fermi surfaces with Ni-dopings of Xe = 0, 0.1, 0.3. The d^^^ dy^, and dxy orbitals for different Fermi surfaces are colored as red, 
green and blue, respectively. (d,e,f) Evolution of low-energy spin excitations for BaFe2-xNi;cAs2 with Xe = 0, 0.1, 0.3. For data in (d) and (f), 
Ei = 450 meV; (e) Ei = 80 meV all with c-axis along incident beam direction. 



In Figure 2 of the main text, we have shown that the widths of spin excitations change with increasing energy above Tc 
for Bao.67Ko.33Fe2As2. The resonance width also changes across Tc (Figure 4 of the main text). To further illustrate the line 
shape change between the normal and superconducting states, we need to determine the width ratio of spin excitations along 
the two high symmetry directions as a function of increasing energy above and below Tc. We plot the fitting of full width at 



half maximum (FWHM) of spin excitations along H and K directions in the normal state in SFig. 11 (a). It is a hour-glass like 
dispersion along [H,0] direction, but a linear dispersion along [1,^] direction. On cooling below T^ the dispersion changes 



to a shape of flower vase in SFig. 11 (b). we plot in SFig. 11 (c) the energy dependence of the ratio (K - H)/(K -h H)fwhm 
from 5 to 36 meV. In the normal state at T = 45 K, spin excitations at energies below the resonance have an oval shape with 
an elongated H direction, thus giving negative anisotropy ratio for energies below about 15 meV. On moving to the resonance 
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SFig 6: Comparison of spin excitations in absolute units for BaFe2-jcNi;,As2 with Xe = 0.1, 0.3. Evolution of low-energy spin excitations 
for BaFe2-^Ni^As2 with = 0.1, 0.3. (a) Data was obtained with Ei = 30 meV for = 0.1 and Ei = 80 meV for = 0.3. (b) Ei = 80 meV 
for X, = 0.1 and Et = 250 meV for = 0.3. (c) Ei = 80 meV for = 0.1 and Ei = 450 meV for = 0.3. (d) Ei = 80 meV for = 0.1 and 
Ei = 250 meV for Xe = 0.3. (e,f,g,h) The corresponding cuts along the [1, K] direction for these two samples. 



energy dXE = 15 meV, the scattering become isotropic and the isotropic scattering persist up to 30 meV. For energies above 30 
meV, transverse elongated scattering take over, giving positive anisotropy ratio. In our previous work |4|, we reported that spin 
excitations integrated from 10 to 18 meV display an longitudinally elongated oval shape in the normal state. This is consistent 
with present work, which has much better energy resolution and better statistics for measured spin excitations. On cooling to 
the superconducting state, scattering are essentially isotropic below the resonance, and change to the maximum anisotropy at the 
resonance energy. For energies above the resonance, the scattering change back to transverse elongated spin excitations above 
35 meV. 

In the main text, the integrated dynamic (local) susceptibilities for hole-doped Bao.67Ko.33Fe2As2, electron-doped 
BaFe1.9Nio.1As2, and BaFe1.3Nio.3As2 are shown in Fig. Ih. To compare spin excitations of Bao.67Ko.33Fe2As2 with 
BaFe2-xNi;cAs2 at Xe = 0,0.1, we show in SFig. [12] two dimensional constant-energy images at different energies in recip- 
rocal space for these materials. The magnetic scattering intensities for the three compounds at the same energy are normalized 
to absolute units in the same color scale. The measurements for Bao.67Ko.33Fe2As2 and BaFe1.9Nio.1As2 were carried out on the 
MERLIN time-of-flight (TOF) chopper spectrometer at the Rutherford- Appleton Laboratory (RAL), UK. The BaFe2As2 mea- 
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SFig 7: Spin excitations in absolute units for BaFei 7Nio 3As2. (a-1 
zone boundary. Data taken on MAPS at 5 K with Ei - 250 meV for (b,( 
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Evolution of spin excitations for BaFe1.7Nio.3As2 all the way to the 
and Ei = 450 meV for (a,c,e,f,g,h). 



surements were carried out at MAPS, RAL. The incident beams were set to be parallel to the c-axis of the sample in all three 
experiments. Although spin excitation energies are coupled to momentum transfers along the c-axis in this scattering geometry, 
we note that MERLIN and MAPS experiments have almost identical L for spin excitation energy if the incident beam energy 
is the same. Therefore, the raw data with the same Ei for Bao.67Ko.33Ee2As2, BaEei.9Nio.i As2, and BaEe2As2 can be compared 
directly. Eor energy transfer of E = 15 + 1 meV, the incident beam energy is Ei = 35 meV for Bao.67Ko.33Ee2As2, 30 meV 
for BaEe1.9Nio.1As2, and 80 meV for BaEe2As2. To compare the data quantitatively, the scattering intensity in 12 a,e,i) have 
been corrected by the magnetic form factor. Spin excitations for Bao.67Ko.33Ee2As2 have the strongest intensity at the resonance 
energy (E = 15 meV) and are rotated 90° in reciprocal space from the line-shape for BaEe1.9Nio.1As2. Spin waves in BaEe2As2 
has a clear anisotropy spin gap below 15 meV. Eor spin excitations at £^ = 50 ± 5 meV in SEig. 12 the incident beam energy was 
Ei = 80 meV. Here, spin excitations for Bao.67Ko.33Ee2As2 are transversely elongated but broader in reciprocal space compared 

[T2I). The scattering intensity is larger than that of BaEe1.9Nio.1As2, 



to those of BaEe1.9Nio.1As2 and BaEe2As2 (SEig. [T2p,[T2f 
but smaller than spin waves of BaEe2As2. However, the integrated dynamic susceptibility of Bao.67Ko.33Ee2As2 is similar with 
that of BaEe2As2, but larger than that of BaEe1.9Nio.1As2. At energies above E > 100 meV, spin excitations of Bao.67Ko.33Ee2As2 
are clearly weaker in intensity than that of BaEe1.9Nio.1As2 and BaEe2As2. The constant-energy images have been subtracted by 
a radial background for Bao.67Ko.33Ee2As2 and a constant background for BaEei.9Nio.i As2 and BaEe2As2. 

To quantitatively determine the dispersions of spin excitations in Bao.67Ko.33Ee2As2, BaEe1.9Nio.1As2, and BaEe2As2, we cut 
through [//, 0] and [1,^] directions of the two dimensional scattering images in SEig. 



ig SEig. [T3^-[T3J: and[T3p-[T3^ show 
constant-energy cuts at energies of E = 15 ± 1, 45 ± 5, 70 ± 5 meV along the [H, 0] and [1, K] directions for Bao.67Ko.33Ee2As2, 
BaEe1.9Nio.1As2, and BaEe2As2, respectively. SEig. 13 1 13 i cuts have been corrected the effect of Bose population factor and 
magnetic form factor, since they have different incident beam energies and T = 45 K data. The cuts along the [1, K] direction at 
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SFig 8: The dispersion of spin excitations for BaFe2-jcNi;tAs2 along [1, ^] direction, (a, c) The dispersion cuts of BaFe2As2 (x = 0) with 
Ei = 250 meV and Ei = 450 meV along [1, ^] direction. The data is from MAPS, (b, d) Identical dispersion cuts of BaFe1.7Nio.3As2 (Xe = 0.3) 
at MAPS. 




SFig 9: Spin excitations cuts in absolute units for BaFe2-xNi;tAs2. (a-f) Evolution of spin excitations along the [1,^] direction for 
BaFe2-xNi;tAs2 at Xe = 0, 0.1, 0.3. (g-i) Similar cuts along the [H, 0] direction. 



E = 135 ± 10, 155 ± 10 and 195 ± 10 meV reveal the weaker susceptibility at high energies for Bao.67Ko.33Fe2As2. The dashed 



blue lines in SFig. [T3g,[T3]i indicate that spin excitations of Bao.67Ko.33Fe2As2 disperse more rapidly and have a softened band 
top. 

Spin excitations disappear at the zone boundary [1, 1] in Bao.67Ko.33Fe2As2, BaFe1.3Nio.3As2, and BaFe2As2. The band top is 
governed by the effective magnetic exchange couplings / (Jia, Jib and J2), as defined in BaFe2As2 |2|. To estimate the change 
of / for hole-doped Bao.67Ko.33Fe2As2, we calculate the energy cut at [1, 1] by exploring the Heisenberg Hamiltonian of parent 
compound. It turns out that /i^. Jib and J2 have comparable eff'ect on the band top. Based on the dispersion of Bao.67Ko.33Fe2As2, 
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SFig 10: Spin excitations cuts in absolute units for BaFe2-jcNi;,As2. (a,b) Constant- 2 cuts of spin excitations for BaFQ2-x^h^^2 at 
Xe = 0, 0.1, 0.3. (c) Coherence lengths of spin excitations for these three samples. 
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SFig 11: Dispersion and anisotropy ratio of the low energy spin excitations for Bao 67Ko 33Fe2As2 above and below T^. (a)Normal and 
(b) superconducting states spin excitation dispersion in full- width at half maximum (FHWM) along the [H,0] and [1,^] directions, (c) The 
filled red circles show the anisotropy ratio at T = 45 K, and the filled blue squares are the same ratio at 9 K. The vertical error bars indicate 
the statistical errors of one standard deviation. The solid lines are a guide to the eye. 



the eflTective magnetic exchange / is found to be about 10 % smaller for Bao.67Ko.33Fe2As2 compared with that of BaFe2As2 
(SFig. [14]). For comparison, if we assume the band top for KFe2As2 is around E = 25 meV, the effective magnetic exchange 
should be about 90% smaller for KFe2As2. Of course, we know this is not an accurate estimation since spin excitations in 
KFe2As2 are incommensurate and have an inverse dispersion. In any case, given the zone boundary energy of £^ ^ 25 meV, the 
effective magnetic exchange couplings in KFe2As2 must be much smaller than that of BaFe2As2. SFigure 15 shows additional 
data for KFe2As2 that clearly reveal a dramatic reduction in magnetic scattering above E = 20 meV. 



RANDOM PHASE APPROXIMATION (RPA) AND DFT+DMFT CALCULATIONS 

As discussed in the main text, RPA calculation of the wave vector dependence of spin excitations in hole-doped 
Bao.67Ko.33Fe2As2 is in clear disagreement with experiments (Fig. 3 of the main text). Although a combined DFT and DMFT 
approach l \ still does not agree in detail with the data (Fig. 3), it captures the trend of spectral weight transfer away from 
Qaf = (1, 0) on increasing the energy and forming a pocket centered at 2 = (1, 1). The solid red, black, and blue lines in SFigure 
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SFig 12: A comparison of constant-energy images of spin excitations for Bao 67K0 33Fe2As2, BaFei qMq 1 As2, and BaFe2As2 as a function 
of increasing energy at low- temperature. The color bars represent the vanadium normahzed absolute spin excitation intensity in the units of 
mbarn/sr/meV/f.u. (a) £ = 15 + 1, (b) £ = 50 + 5, (c) £ = 155 + 10, and (d) £ = 195 + 10 meV are for Bao.67Ko.33Fe2As2 at 9 K. (e-h) and 
(i-1) are identical images for BaFei.gNio.i As2 at 5 K and BaFe2As2 at 7 K, respectively. 



16] show calculated local susceptibility in absolute units based on a combined DFT and DMFT approach for Bao.67Ko.33Fe2As2, 



BaFe2As2, and BaFe1.9Nio.1As2 |6|, respectively. This theoretical method predicts that electron doping to BaFe2As2 does not 
affect the spin susceptibility at high energy {E > 150 meV), while spin excitations in the hole doped compound beyond 100 
meV are suppressed by shifting the spectral weight to lower energies. This is in qualitative agreement with our absolute intensity 
measurements (Fig. Ih). The reduction of the high energy spin spectral weight and its transfer to low energy with hole doping, 
but not with electron doping, is not naturally explained by the band theory, and requires models which incorporate both the 
itinerant quasiparticles and the local moment physics. The hole doping makes electronic state more correlated, as local moment 
formation is strongest in the half-filled shell, and mass enhancement larger thereby reducing the electronic energy scale in the 
problem. 

Our theoretical DFT+DMFT method for computing the magnetic excitation spectrum employs the abinitio full potential 
implementation of the method, as detailed in O. The DFT part is based on the density functional theory (DFT) code of Wien2k 
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SFig 13: Constant-energy cuts of spin excitations as a function of increasing energy, (a-c) are the cuts along the [H,0] direction and 
(d-i) are along the [1,^] direction for Bao.67Ko.33Fe2As2 at 9 K (red filled circles), 45 K (half filled purple squares), BaFe1.9Nio.1As2 at 5 K 
(solid blue line) and BaFe2As2 at 7 K (solid black line), (sl) E = 15 + 1 meV with Ei = 35 meV for Bao.67Ko.33Fe2As2 and Ei = 30 for 
BaFe1.9Nio.1As2, (b) £ = 45 + 5 meV with Ei=80meV, (c) £ = 70 + 5 meV with Ei = 250 meV are cuts along [H,0] direction, (d-f) are the 
same constant-energy cuts along [1, ^] direction, (g) E = 135 + 10, (h) £^ = 155 + 10, (i) £^ = 195 + 10 meV are from Ei = 450 meV data. 
The solid lines are Gaussian fits to the data. (a,d) have been corrected the Bose population factor and magnetic form factor. 
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SFig 14: The effect of magnetic exchange couplings / (/i^, Jit and J2) on the band top of spin excitations. The black line is energy cut at 
(0.8 <H <1.2, 0.8 <H <1.2) r.l.u for BaFe2As2 in the Heisenberg spin wave model |2|. The red line is for Bao.67Ko.33Fe2As2 with 10 % soften 
band top. The blue line is a similar estimation for KFe2As2 assuming zone boundary is around E = 25 meV. 



191 . The DMFT method requires solution of the generalized quantum impurity problem, which is here solved by the numerically 
exact continuous-time quantum Monte Carlo method |[TOl[TT]|. The Coulomb interaction matrix for electrons on iron atom was 
determined by the self-consistent GW method in Ref. lfT2ll . giving U = 5 and / = 0.8 eV for the local basis functions within 
the all electron approach employed in our DFT-fDMFT method. The dynamical magnetic susceptibility ;^''( 2, E) is computed 
from the ab initio perspective by solving the Bethe-Salpeter equation, which involves the fully interacting one particle Greens 
function computed by DFT-hDMFT, and the two particle vertex, also computed within the same method (for details see Ref. 
fTl). We computed the two-particle irreducible vertex functions of the DMFT impurity model, which coincides with the local 
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SFig 15: The disappearance of sin excitations above 25 meV for KFe2As2. Using Ei - 80 meV, we can still see clear incommensurate spin 
excitations at (a) E - 17.5 + 2.5 meV. However, for energy transfers dXE - 22.5 + 2.5 mV, one can no longer see spin excitations, indicating a 
strong suppression of the spin excitation spectral weight for E > 25 meV. 



two-particle irreducible vertex within DFT+DMFT method. The latter is assumed to be local in the same basis in which the 
DMFT self-energy is local, here implemented by projection to the muffin-tin sphere. 
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SFig 16: RPA and LDA+DMFT calculated local susceptibility for different iron pnictides. RPA and LDA+DMFT calculations of x"{oj) 
in absolute units for Bao.67Ko.33Fe2As2 comparing with earlier results for BaFe2As2 and BaFei.gNio.i As2 | 6 |. 



MAGNETIC EXCHANGE ENERGY AND SUPERCONDUCTING CONDENSATION ENERGY FOR BAq 67K0 33FE2AS2 

In a neutron scattering experiment, we measure scattering function S(q,E = ficS) which is related to the imaginary part of 
the dynamic susceptibility VmS{q, cS) = [I + n(aj, T)]x'\q, co), where [1 -h n(<jj, T)] is be Rose population factor. The magnetic 
exchange coupling and the imaginary part of spin susceptibility are related via the formula |[T3l : 



{Si'S j) = [ -^f-r r dco[\ + n(oj, T)]x\q, oj)cos[q • (/- /)], (1) 

^8 J^B J ^^^^ J 
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where g = 2 is the Lande g-factor. The magnetic exchange energy can be written as 



<i,J> 

3 ^ 



i 

+ ^ + n{(jj, T)]y\q, oj)cos(qy) 

i 

+ ^ /2[1 + n(oj, T)]x'\q, oj)[cos(qx + qy) + cosiq^ - qy)]}. 



(2) 



Here we have assumed an anisotropic model for the effective magnetic exchange coupHng |2|, different from the case of copper 
oxide superconductors (13]. Jia is the effective magnetic coupHng strength between two nearest sites along the a direction, 
while Jii) is that along the b direction, and J2 is the coupling between the next nearest neighbor sites. Hence we are able to 
obtain the change in magnetic exchange energy between the superconducting and normal states by the experimental data of 
X'\q, ^) in both states. Strictly speaking, we want to estimate the zero temperature difference of the magnetic exchange energy 
between the normal and the superconducting states, and use the outcome to compare with the superconducting condensation 
energy |[T3ll . Unfortunately, we do not have direct information on the normal state ^) at zero temperature. Nevertheless, 
since our neutron scattering measurements at low-energies showed that the x'\q^ ^) are very similar below and above Tc near 
the AF wave vector Qaf = (1, 0, 1) and only a very shallow spin gap at 2 = (1, 0, 0) (see Figs. If and Ih in |4 |), we assume that 
there are negligible changes in;^'^^. ^) above and below Tc at zero temperature for energies below 5 meV. For spin excitation 
energies above 6 meV, Bose population factors between 7 K and 45 K are negligibly small. In previous work on optimally doped 
YBa2Cu3 06.95 superconductor, we have assumed that spin excitations in the normal state at zero temperature are negligibly small 
and thus do not contribute to the exchange energy llT4l . 

The directly measured quantity is the scattering differential cross section 



d cr k; 



dD.dE k 



2(yr )^ ^ 

^l^\F(Q)f[l^n(oj,T)]y\q,oj), 



(3) 



where ki and kf are the magnitudes of initial and final neutron momentum and F(Q) is the Fe magnetic form factor, and 
(yr,)2 = 0.2905 barn-sr-i. 



The quantity 22 x'^iq^ ^) in both superconducting and normal states can be fitted by a Gaussian As{n)e 
resonance wave vector (1,0) and by cutting the raw data. The outcome is summarized in the table [] 



for 



E 


0-x,s 




As 


0-x,n 




An 


5 


0.050 


0.060 


6.017 


0.115 


0.076 


6.571 


1 


0.059 


0.066 


7.318 


0.109 


0.080 


5.929 


9 


0.077 


0.077 


9.789 


0.154 


0.100 


5.515 


11 


0.092 


0.083 


12.001 


0.160 


0.107 


5.745 


13 


0.121 


0.097 


14.674 


0.145 


0.124 


5.102 


15 


0.153 


0.106 


16.792 


0.125 


0.130 


6.535 


17 


0.167 


0.118 


12.141 


0.152 


0.157 


5.546 


20 


0.173 


0.136 


8.856 


0.131 


0.128 


5.555 


24 


0.165 


0.155 


6.802 


0.134 


0.146 


5.411 


28 


0.182 


0.177 


3.393 


0.161 


0.165 


3.393 



where the unit of E is meV and that of As(n) is mbam-meV"^ -sr"^ -Fe'^ For the case below 5 meV, we assume that A„ decreases 
to zero linearly with energy and As = A„ (see Fig. Ih in Ref. O), while the cr's keep the values at 5 meV. The assumption is 
shown in SFig. [] where the resonance is seen at E = 15 meV. 
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Because the condensation energy is only defined at zero temperature, we take T = in equation (1) and the integral gives: 

{Si ■ 5,+,), - {Si ■ Si^,}„ = -0.0039, 
{^i ■ ^i^y)s - {^i ■ Si^y)n = 0.0043, 

J] «^,- • §i^l)s - {Si ■ Si^l)n) = -0.0073 

l^x±y 

The magnetic exchange couphng constants in an anisotropic model are estimated to be 

JiaS = 53.3 meV, 

JitS = -8.3 meV, (5) 
J2S = 12.2 meV, 

which are 10% smaller than that of BaFe2As2 | 2 | and we estimate S to be close to \ ESIITtI. Hence the exchange energy change 
is 

^Eex = -0.66 meV/Fe. (6) 
The condensation energy Uc for optimally doped Bao.68Ko.32Fe2As2 can be calculated to be 

Uc = -17.3 J/mol 

.^o lev 1 
= -17.3 



1.6 X 10-19 6.02 X 1023f.u. 

_ lev 1 

1.6 X 10-19 2 X 6.02 X 1023 Fe 
= -0.09 meV/Fe 

from the specific heat data of Ref. (15]. Therefore, we have the ratio AEex/Uc ~ 7.4, meaning that the change in the magnetic 
exchange energy is sufficient to account for the superconducting condensation energy in Bao.68Ko.32Fe2As2. We note that a 
similar calculation for heavy Fermion superconductor CeCu2Si2 also reveals that the change in magnetic exchange energy is 
sufficient to account for the superconducting condensation energy [TSl . 

We thank M. S. Liu and L. W. Harriger for providing the BaFe1.9Nio.1As2 and BaFe2As2 raw data, respectively. 
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